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Abstract. We generalize the supersymmetry method in Random Matrix Theory to 
arbitrary rotation invariant ensembles. Our exact approach further extends a previous 
contribution in which we constructed a supersymmetric representation for the class 
of norm-dependent Random Matrix Ensembles. Here, we derive a supersymmetric 
formulation under very general circumstances. A projector is identified that provides 
the mapping of the probability density from ordinary to superspace. Furthermore, 
it is demonstrated that setting up the theory in Fourier superspace has considerable 
advantages. General and exact expressions for the correlation functions are given. 
We also show how the use of hyperbolic symmetry can be circumvented in the 
present context in which the non-linear a model is not used. We construct exact 
supersymmetric integral representations of the correlation functions for arbitrary 
positions of the imaginary increments in the Green functions. 



PACS numbers: 05.45.Mt, 05.30.-d, 02.30.Px 
1. Introduction 

The supersymmetry method is nowadays indispensable for the discussion of various 
advanced topics in the theory of disordered systems [H and it became equally 
important in numerous random matrix approaches to complex systems in general |31 
m El ini- Random Matrix Theory (RMT) as originally formulated in ordinary space 
does not rely on Gaussian probability densities. It is only important that the Random 
Matrix Ensembles are invariant under basis rotations. Gaussian probability densities 
are highly convenient in calculations, but other probability densities are also possible, 
and some of those were already considered in the early days of RMT [Zj. On the 
other hand, the supersymmetric formulations were constructed for Gaussian probability 
densities [H El E] by means of a Hubbard-Stratonovich transformation. Thus, the 
question arises naturally whether the Hubbard-Stratonovich transformation restricts 
the use of supersymmetry to the Gaussian form of the probability densities. We address 
this problem in the present contribution. We will show that the supersymmetry method 
is not at all restricted in this way, and we will derive supersymmetric formulations of 
RMT for arbitrary rotation invariant Random Matrix Ensembles. 
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We focus on conceptual and structural issues. In particular, we are not aiming at 
asymptotic results in the inverse level number as following from the supersymmetric non- 
linear a model [HEIIS]- This latter approach was used in Ref. P to show universality for 
infinite level-number in the case of non-Gaussian probability densities. Here, however, 
our goal is different: we address the full problem to achieve exact, i.e. non-asymptotic 
results. In a previous study ^U], we presented supersymmetric representations for norm- 
dependent ensembles, where the probability densities are functions of the traced squared 
random matrices only. Although a series of interesting insights are revealed already 
in this case, the derivation can be done without actually employing deep features of 
supersymmetry. This is not so in the present contribution which aims at a general 
construction. The methods needed are very different from the ones of Ref. Here, 
we have to explore the algebraic structure of superspace. 

One can also motivate the present investigation from the viewpoint of applications. 
We refer the interested reader to the contribution ^| and the literature quoted therein. 
Our goal to perform a conceptual study does not prevent us from giving general 
expressions for the correlation functions, but we refrain from looking too much into 
applications and defer this aspect to future work. 

It will not be surprising for those who already have expertise in supersymmetry 
that a generalization as outlined above requires an analysis of convergence properties 
and thus leads inevitably to the issue of what kind of symmetries the theory in 
superspace should have. It was argued in Ref. that hyperbohc symmetry, i.e. groups 
comprising compact and non-compact degrees of freedom, are necessary if one is to 
set up a non-linear a model in ordinary space. This line of reasoning carries over to 
superspace [HE], see also the recent review in Ref. jH]. We justify a procedure for how 
to avoid hyperbolic symmetry in the framework of our supersymmetric models. The 
necessity to introduce hyperbolic symmetry is exclusively rooted in the non-linear a 
model, not in supersymmetry as such. If i.e. non-asymptotic results, 

compact supergroups suffice. 

For various reasons, including some related to convergence questions, we find it 
advantageous to map the theory onto Fourier superspace. Moreover, we restrict ourselves 
to unitary Random Matrix Ensembles throughout the whole study. 

The paper is organized as follows. Having posed the problem in Section |2l we 
generalize the Hubbard-Stratonovich transformation in Section 01 In Section HI we 
derive the supersymmetric formulation in Fourier superspace. The correlation functions 
are expressed as eigenvalue integrals in Sectional Summary and conclusions are given 
in Section ini 

2. Posing the Problem 

In Section 12.11 the two relevant kinds of A;-point correlation functions are defined and 
the relation to the generating functions is given. Thereby we also introduce our notation 
and conventions. We clarify what we mean by arbitrary rotation invariant ensembles in 
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Section 12.21 In Section 12.31 we show how different types of correlation functions can be 
related to each other by proper Fourier transforms. 

2.1. Correlation and Generating Functions 

The Random Matrix Ensemble builds upon N x N Hermitean matrices H, having 
altogether iV^ independent matrix elements. A normalized probability density P{H) 
assigns a statistical weight to the elements of the matrices H. As the Hermitean matrices 
are diagonalized by unitary matrices in SU(A^), the probability density P{H) is said to 
define a Unitary Random Matrix Ensemble. We are interested in the fc-point correlation 
functions 

k 

Rk{x,,...,Xk)= J d[H] P{H) n tr 5 {x^ - H) , (1) 

p=i 

depending on the k energies xi,...,Xk- The 6 functions are the imaginary parts of 
the matrix Green's functions, ^inS {xp — H) = Im {xp±ie — H)~^. Here, ie is an 
imaginary increment and the limit £ ^ is suppressed in the notation. In the 
supersymmetric construction to follow, it is convenient to consider the more general 
correlation functions which also include the real parts of the Green's functions. They 
are, apart from an irrelevant overall sign, given by 

Mxu ...,x,) = ^f d[H] P{H) n tr . (2) 

One often wants to put the imaginary increments on different sides of the real axis. 
The quantities Lp which are either +1 or —1 determine the side of the real axis where 
the imaginary increment is placed. The correlation function can always be expressed as 
derivatives of a generating function Zk{x + J) such that 



Rkixi, . . . , Xk) = -, n ; Zuix + J) 



(3) 



where 



We introduced source variables Jp, p = 1, . . . , as well as the diagonal matrices 
X = diag (xi, xi, . . . , Xk, Xk) and J = diag (— Ji, + Ji, . . . , —Jk, +Jk)- In the sequel, we 
use the short hand notations Xp = Xp — iLpS and x^ = diag (xf , xf , . . . , x^, x^). The 
product of the differentials of all independent matrix elements is the volume element 
d[H]. We use the notation and the conventions of Refs. 11^1 IT^. The normalization 
Zk{x) = 1 at J = follows immediately from the definition (@)). 

We wish to study whether the generating function can be represented as an integral 
of the form 

Zfc(x + J) = j d[a] g(a)detg [a - x^ - j) , (5) 
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where a is a 2k x 2k supermatrix with Hermitean or related symmetries, and where 
detg denotes the superdeterminant. If such a representation can be shown to exist, the 
question arises whether the probabihty density Q{ct) in superspace can be obtained in 
a unique way from the probabihty density P{H) in ordinary space. 

2.2. Rotation Invariant Probability Densities 

For the important class of norm-dependent ensembles, i.e. ensembles defined by a 
probability density depending exclusively on trif^, such a unique construction is indeed 
possible and was performed in Ref. fOj- Here, we tackle the problem of arbitrary 
rotation invariant probability densities P{II). We recall that a probability density must 
be normalizable and positively semi-definite. The term "arbitrary" has to be understood 
as excluding those functions P{II) which would lead to a divergent integral By 
"rotation invariant" we mean that the probability density has the property 

P{H) = PiUoHUl) = P{E) , (6) 

where Uq is any fixed matrix in SU(A^) and where E = diag {Ei, . . . , Ep^) is the diagonal 
matrix of the eigenvalues of H. Although it is obvious, we underline that this includes 
invariance under permutations of the vectors defining the basis in which H is written 
down and also invariance under permutations of the eigenvalues. Hence, the probability 
density P{II) should depend only on matrix invariants, such as ti H"^ where m is 
real and positive. Anticipating the later discussion, we already now mention that this 
requirement is a most natural one in view of the general character of the supersymmetry 
method. The strength of this method is rooted in the drastic reduction of degrees of 
freedom, i.e. of the number of integration variables, when an integral over the N x N 
matrix H is identically rewritten as an integral over the 2k x 2k matrix a. Thus, 
supersymmetry removes a certain redundancy. The rotation invariance requirement 
implies precisely this redundancy which the supersymmetry method needs. We will 
show that this holds for arbitrary rotation invariant probability densities P{II). 

2.3. Mutual Relations between the Different Correlation Functions 

We wish to address the correlation functions and for finite level number A^, we 
are not aiming at an asymptotic discussion. If a saddlepoint approximation leading 
to a non-linear a model as in Refs. PJ El E] is the method of choice to study a 
certain physics problem, one performs precisely such an asymptotic expansion in 1/A^. 
This is not what we are going to do in the present contribution. Admittedly, our 
goal to address the problem exactly for finite A^ renders our task mathematically 
demanding, because we have to solve certain group integrals. One the other hand, 
luckily and at first sight paradoxically, this goal allows us to circumvent the introduction 
of hyperbolic symmetry, which is a deeply rooted, non-trivial feature of the non-linear 
a model pi El El E]- Hyperbolic symmetry means that the ensuing supersymmetric 
representation of the random matrix model must involve non-compact groups to make 
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the integrals convergent. This is inevitable if the imaginary increments of the energies 
lie on different sides of the real axis. However, if they lie on the same side, no hyperbolic 
symmetry occurs and all groups are compact. This facilitates the supersymmetric 
treatment tremendously. 

We now argue that the correlation functions (P) of the imaginary parts can be 
recovered from the more general correlation functions (j21) that are suited for the 
supersymmetric treatment, even if all imaginary increments lie on the same side of 
the real axis. We choose Lp = +1 for all p = 1, . . . ,k. Upon Fourier transforming the 
correlation functions (0), we obtain the /c-point correlations in the domain of the times 
tp-i P 1, . . . , fc, 

^ +00 +00 

ffe(ti, . . . , tk) = ^=fc / dxi exp (itiXi) ■■■J dxk exp (it^Xk) 

— oo 

Rk{xi, . . . , Xk) 



271 



= it2y n exp i-etp) rfc(ti, ...,tk) (7) 
p=i 

1 ^ 
Tkiti, ...,tk) = I d[H] P{H) n trexp {iHtp) . (8) 

v27r •' p=i 



with 



Importantly, this latter fc-point correlation function rk{ti, . . . ,tk) in time domain is 
precisely the Fourier transform of the correlation function (^. It is well-defined on the 
entire real axes of all its arguments tp. The inverse transform yields 

^ +00 +00 

Rkixi, ...,Xk)= / dti exp {-ixiti) ■ ■ ■ / dt^exp {-ixktk) 



-OO 



rk{ti, ...,tk) . (9) 

Based on this observation, we will pursue the following strategy in later Sections of this 
contribution: We perform exact manipulations of the correlations Rk{xi, . . . ,Xk) with 
Lp = +1, p = 1, . . . , fc, or of their generating functions, respectively. Having obtained 
the appropriate supersymmetric representation, we Fourier transform it into the time 
domain and find ffc(ti, . . . ,tk). In this expression, we then identify the supersymmetric 
representation of the correlation functions rk{ti, . . . ,tk). Upon backtransforming 
we arrive at the desired supersymmetric representation for the correlation functions 
Rk{xi, . . .,Xk). 

We can even extend the line of arguing. Once we have found the supersymmetric 
representation of rk(ti, . . . ,tk), we can construct the one of Rk{xi, . . . ,Xk) for any 
arbitrary choice of the quantities Lp = ±1 by calculating 

^ +00 +00 

-Rfc(xi, ...,Xk) = / dti exp {-ixiti) ■ ■ ■ / dtkexp {-ixktk) 



27T 

k 

n (^''^2) Q{Lptp) exp i-eLptp) r^it^, . . . , 4) • (10) 



oo 

k 



p=l 
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In this manner, we will obtain supersymmetric integral representations for all correlation 
functions (0) where the imaginary increments lie independently of each other on either 
side of the real axes, without introducing hyperbolic symmetry. 

3. Generalized Hubbard— Stratonovich Transformation 

To carry out the program outlined in Section |21 we have to generalize the procedure 
referred to as Hubbard-Stratonovich transformation accordingly. In Section 13. H 
we Fourier transform the probability density. An algebraic duality between matrix 
structures is uncovered in Section ing and explored further in SectionESl where spectral 
decompositions of the matrices involved are performed. Although our main interest 
are the correlation functions where all Lp are equal, we make these latter steps for an 
arbitrary metric L. We do so, because we find it worthwhile to document how natural the 
duality is even for a general metric. Moreover, it allows us to clearly identify the point 
where a general metric would require a much involved discussion of hyperbolic symmetry 
— which we then avoid by setting Lp = +1 for all p = 1, . . . ,k. In Section \liA\ we 
construct the probability density in superspace. We derive a generalized transformation 
formula and the corresponding generating function in Sections 13 .51 and IH .fi| respectively. 
In Section EIZl the norm-dependent ensembles are discussed as a simple example. 

3.1. Fourier Transform of the Probability Density 

The determinants in the generating function (j3)) are written as Gaussian integrals, 
those in the denominator as integrals over k vectors Zp, p = 1, . . . , k with complex 
commuting elements each, and those in the numerator over k vectors (p, p = 1, . . . ,k 
with A^ complex anticommuting elements each. Again omitting irrelevant phase factors, 
we have 



where d[zp] and d[Cp] denote the products of the independent differentials. To ensure 
convergence of the integrals over the commuting variables, the quantities Lp are inserted 
in front of the bilinear forms in the exponent. This is not needed in the integrals over the 
anticommuting variables because they are always convergent. We order the quantities 
Lp in the metric tensor L = diag (Li, 1, . . . , L^, 1). Using the identities 




Jl / d[zp]exp [iLpzl{H 






zlHzp = tr Hzpzl and QHCp = - tr HCpCp , 
the average over H in Eq. (fTT|) can be written as the Fourier transform 



(12) 




(13) 
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of the probability density. The Fourier variable is the matrix 



k k 



(14) 



p=i p=i 



The function ^{K) is referred to as characteristic function. The definition (fTSj) of the 
Fourier transform is the one mostly used in the statistics literature. It guarantees that 
$(0) = 1, directly reflecting the normalization of P{H). The definition of the Fourier 
transform in Section 12.31 follows the "symmetric convention" in which the same factor 



Up to now, all steps were exactly as in the case of a Gaussian probability density 
P{H). In the Gaussian case, one can now do the integral ()13|) explicitly and one obtains a 
Gaussian form for the characteristic function (^{K). Here we consider a general rotation 
invariant P{H). Of course, we must assume that the Fourier transform exists, i.e. that 
P{H) is absolutely integrable or, even better, that it is a Schwartz function. Absolute 
integrability is guaranteed by the fact that P{H) is a probability density, implying that 
it is positively semi-definite and normalized. However, we also must assume that the 
integrals over the vectors Zp converge after doing the Fourier transform. The integrals 
over the vectors (p can never cause convergence problems. In the Gaussian case, all 
those convergence issues have been carefully discussed in Ref. fl], a recent review is 
given in Ref. P|. In the general rotation invariant case, we have no other choice than 
to implicitly exclude those probability density P{H) which would cause convergence 
problems, assuming that all integrals in the sequel converge. We will come back to this 
point later. 

It is easy to see that the rotation invariance of P{H) also implies the rotation 
invariance of ^{K). The matrix K is Hermitean, K'' = K. This is so for all choices 
Lp = ±1 of the metric elements. As the entries of K are commuting variables, we may 
conclude that K can be diagonalized. 



where V is in SU(A^) and where Y = diag {Yi, . . . , Yn) is the diagonal matrix containing 
the eigenvalues Yn, n = 1,...,N of K. The rotation invariance of P{H) and the 
invariance of the measure d[H] allows one to absorb V such that the characteristic 
function $(i^) depends only on Y, 



In other words, ^{K) is a rotation invariant function, too. 
3.2. Underlying Algebraic Duality 

The merit of the Hubbard-Stratonovich transformation in the supersymmetry method 
is the drastic reduction in the number of degrees of freedom. This is rooted in a duality 
between matrices in ordinary and superspace. We uncover this duality and cast it into 




K = VYV^ , 



(15) 




(16) 



Arbitrary Rotation Invariant Random Matrix Ensembles and Supersymmetry 



8 



a form which allows a straightforward generalization of the previous discussion for a 
Gaussian probability density. We define the N x 2k rectangular matrix 

A=[z,---ZkCi--<k] . (17) 

Although it contains commuting and anticommuting entries, A is not a supermatrix of 
the type commonly appearing in the framework of the supersymmetry method Pfl 
Nevertheless, this matrix will play a crucial role in the following. Its Hermitean 
conjugate is the 2k x N rectangular matrix 



At 



The inclusion of the minus signs is necessary to be consistent with the conventions in 
Refs. UmailSllIl. It ensures that we have {A^Y = A. 

We notice that the boson-fermion block notation |3] is used in the definition p7|l . 
which differs from the pq block notation [3| employed when defining the supermatrices L, 
X and J, as well as implicitly a in Section |2l In the boson-fermion block notation, first 
all commuting and then all anticommuting variables (or vice versa) are collected in a 
supervector. Hence, the supermatrices which linearly transform those vectors consist of 
rectangular (in the present case kxk) blocks of commuting or anticommuting variables. 
The pq notation is obtained by simply reordering the basis. One collects the commuting 
and anticommuting variables corresponding to each energy index p = 1, . . . ,k, such 
that every supermatrix is written as a k x k ordinary matrix with 2x2 supermatrix 
elements assigned to each index pair {p,q). While the latter notation was handy when 
introducing the generating function in Section |2l it is more convenient for the present 
discussion to use the boson-fermion block notation. In particular, the metric then reads 
L = diag (Li, . . . , L^, 1, . . . , 1). 

The Hermitean N x N matrix K defined in Eq. fll4j) can be written as the matrix 
product 

K = ALA^ = (AL^/^) (L^/^A^) . (19) 

There exists a natural dual matrix to K, found by interchanging the order of the matrices 
in Eq. (fTTHl . It is the 2k x 2k matrix 

B = (L^/^At) {AL'/^) = L'/'A^AL'/' , (20) 
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where 



z[zi 



4^1 
-Clzi 



<l 



-cki 



ZlCk 

ziCk 
-ckk 



-(kCk 



(21) 



-Cfc^i ■ ■ ■ -(Izk -(Id ■ ■ 
While K = ALA^ is an ordinary matrix, A and B = L^^"^ AL^^'^ are supermatrices. 
Moreover, K is Hermitean for all choices of the metric L, i.e. for every combination 
Lp = ±1, but B is in general not Hermitean because some entries of the metric are 
imaginary, L^/^ = i. The supermatrix A'^A, however, is Hermitean. 

Interestingly, the duality between the matrices K and B also implies the equality 
of invariants involving the traces according to 

trK" = trg5™, (22) 

for every non-zero, positive integer m. This generalizes the case of a Gaussian 
probability density where the need to discuss this equality occurs only for m = 2. As 
the equality is not completely trivial due to the presence of anticommuting variables, 
Eq. ()22|1 is proven in [Appendix A 



3.3. Eigenvalues and Eigenvectors of the Dual Matrices 

Our way of formulating the algebraic duality is most helpful for the spectral 
decomposition in ordinary and superspace. We write the eigenvalue equation for the 
matrix K as 



(23) 

with eigenvectors Ki, n = 1, . . . ,N. We will now construct them in such a way 
that they are not identical to those given as the columns Vn of the unitary matrix 
V introduced in the diagonalization ()15|) . For our construction, we employ the 2k 
component supervectors 

Wnll 



Wr. 



(24) 



Wnkl 
Wnl2 

Wnk2 

There are two distinct representations of these supervectors. In the first one, the 
elements Wnpj are commuting if j = 1 and anticommuting if j = 2, in the second 
representation it is the other way around. We make the ansatz 

k k 

Vn = AL^^^Wn = ZpLi^^Wnpl + CpWnp2 (25) 
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for the eigenvectors. It is convenient to multiply the coefficients Wnpj from the right to 
avoid some cumbersome signs if the Wnp2 are anticommuting and appear together with 
the vectors (p. We plug the ansatz (jSHI) into the eigenvalue equation (jSHl) and find 

KVr, = ALA^AL'^/^Wn = AL^^^ Bw^ 

Y^Vn = Y^AL'/^Wn = AL^/^ Y^w^ , (26) 
which yields AL^^'^{Bwn — YnWn) = 0. Hence, we conclude that the eigenvalue equation 

BWn = YnWn (27) 

holds if the eigenvalue equation is valid and if the eigenvectors Vn have the form ()25|l . 
There is a duality: the eigenvalues Yn of K to the eigenvectors Vn in the form ()25|) are 
also eigenvalues of B to the eigenvectors w„. 

The fact that the eigenvectors w„ of the supermatrix B belong to one distinct 
representation as discussed below Eq. implies that there are two types of eigenvalues 
corresponding to these representations. We denote the k eigenvalues associated with the 
ffist representation hj i/pi = Yp, p = 1, . . . , k and the k eigenvalues associated with the 
second one by yp2 = Yk+p, p = 1, . . . ,k, respectively. Moreover, not all eigenvectors 
Vn of K can have the form (^3)) if the vector Wn is required to be eigenvector of B at 
the same time. This is so, because K and B have different dimensions N x N and 
2k X 2k, respectively. In all applications of RMT and supersymmetry, the level number 
N is large, such that we may safely assume > 2k. The matrix B has 2k eigenvalues. 
Thus, the duality uncovered above only makes a statement about 2k out of the 
eigenvalues of K. Importantly, the remaining eigenvalues of K are zero, because K is 
built upon 2k dyadic matrices. Hence, we have 

{ypi for n = p, p = l,...,k 
yp2 for n = p + k, p = 1, . . . ,k , (28) 
for n = 2A; + l,...,A^ 

if > 2k. As K is an ordinary Hermitean matrix, we know that the eigenvectors Vn 
to the zero eigenvalues can be chosen orthogonal with each other and with those to 
the non-zero eigenvalues. We order the non-zero eigenvalues in the 2k x 2k diagonal 
supermatrix 

y = diag {yu, . . . , yki, iyk2, • • • , iyk2) (29) 

in boson-fermion block notation. The definition includes an imaginary unit i coming 
with all eigenvalues yp2. This is done for convenience, the motivation will become clear 
later. 

As the presence of the anticommuting variables requires some care, the line of 
reasoning given above is supplemented with some details in [Appendix B[ including the 
relation between the eigenvectors Vn and Vn- 
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3.4- Probability Density in Superspace 

The characteristic function $(K) of the probabihty density P{H) is according to Eq. |T6|l 
rotation invariant, ^{K) = ^{Y). Furthermore, by virtue of the previous discussion we 
may view it as function of the eigenvalues of the supermatrix B, such that we arrive at 
the chain of equahties 

^{K) = = $(?/) = $(5) . (30) 

This crucial observation identifies $ as an invariant function in two different spaces, in 
ordinary space depending on the N x N matrix K and in superspace, depending on 
the dual 2k x 2k matrix B. It is interesting to notice that, if ^{K) is a function of 
all invariants ti K"^ with m = 1,2,3, .. ., we may conclude from the equality ()22|1 the 
identity 

$(trf^,tri^^trir^...) = <l>(trg5,trg5^trg5^ . . .) , (31) 

implying that the form of $(fC) as function of those invariants fully carries over to 
superspace. Although we have no reason to doubt that this also holds in the presence 
of invariants tvK"^ with non-integer m, we have no proof, because we had to assume 
integer m when deriving the equality (j^^ . Luckily, this is not important in the sequel. 
All what really matters is the general insight expressed by the chain of equalities (jHUj) . 
It includes all invariants, except the determinant det K which is trivially excluded, since 
we know from the previous Section 13.31 that K has zero eigenvalues if > 2k. Thus, 
we do not employ Eq. (pTH) in the sequel, although we will refer to it at one point for 
illustrative purposes. 

We restrict the further discussion to the case that all imaginary increments of the 
energies lie on the same side of the real axis. Hence we choose the metric 

L = +l2k ■ (32) 

This choice implies that the supermatrix B becomes Hermitean, B"^ = B, and the 
symmetry group is \J{k/k), the unitary supergroup in k bosonic and k fermionic 
dimensions. For a general metric, the corresponding relation reads B'^ = LBL. The 
symmetry group is pseudounitary, i.e. the matrices of the defining representation satisfy 
w'^Lw = L. This hyperbolic symmetry involves non-compact degrees of freedom. The 
situation was analyzed in detail in Refs. jTI] and PP for the non-linear a model in 
ordinary and in superspace, respectively. The proper, convergence ensuring integration 
manifolds of the Hubbard-Stratonovich fields, corresponding to the matrices a in the 
present case, was constructed. It seems to us that the Gaussian form of the probability 
densities P{H) in these investigations was somehow important for this construction. 
Here, however, we study arbitrary rotation invariant probability densities P{H). We 
did not succeed in extending the line of reasoning in Refs. PflE] to such general P{H), 
even though we strongly believe that this should be possible. Nevertheless, this does not 
cause a problem in view of what we are aiming at, because we can proceed as outlined 
in Section 12.31 All issues related to convergence can be dealt with much easier if the 
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choice (|32j) is made. That this works fine in the case of a Gaussian P{H) was aheady 
demonstrated in Ref. |12j . 

As is the characteristic function of P{H) in ordinary space, the chain of 

equahties naturally suggests to interpret $(-B) as a characteristic function in 
superspace. To this end, we introduce a probability density Q{(j) depending on a 2A; x 2k 
supermatrix a whose Fourier transform is $(-B). However, there is a subtle point to 
which we have to pay attention. The symmetries of B dictate to a large extent what 
the symmetries of a have to be. As 5 is a Hermitean supermatrix, a ought to be a 
Hermitean supermatrix as well. Nevertheless, a Wick-type-of rotation was applied in 
the case of Gaussian probability densities which provides all elements in the fermion- 
fermion block of a with an imaginary unit z PHH]. This modification is needed to solve a 
convergence problem, too. It makes the expression trg cr^ positive semi-definite, and thus 
the integrals over the Gaussian probability density Q{(j) ~ exp(— trg cr^) convergent. 
As we want to include the Gaussian as a special case in our considerations, we also 
introduce this Wick-type-of rotation in the 2k x 2k supermatrix a. The entries of 
the matrix B can be modified accordingly by multiplying the vectors (p containing the 
anticommuting variables with factors of Vi, if one wishes, but we do not do that here. 
The diagonalization of the matrix a can be written as 

a = usu^ with s = diag (sn, . . . , s^i, zsi2, • • . , zsfc2) , (33) 

where all eigenvalues Spj are real. The Wick-type-of rotation multiplies the eigenvalues 
Sp2 with an imaginary unit. Thus, u is in the unitary supergroup U{k/k), without any 
modification of its matrix elements. We also introduce a. 2k x 2k supermatrix p with 
the same symmetries as a and with the diagonalization 

p = vrv'' with r = diag(rii, . . . ,rfci,iri2, . . . ,irfc2) , (34) 

where v is in the unitary supergroup l]{k/k) as well. Anticipating the definitions of 
the eigenvalue matrices s and r, we introduced the eigenvalue matrix y oi B in the 
form PH) . 

We now define the probability density Qi^a) in superspace through the Fourier 
integral 

J d[a] Q{a) exp (itigar) = $(r) , (35) 

or, as <l>(r) is invariant, we have equivalently 

J d[a] Qia) exp (ttigap) = ^(p) , (36) 

where $(p) is obtained from $(-B) by formally replacing B with p. We recall that the 
matrix elements of B in the fermion-fermion block are the scalar products CpCg ^iid thus 
nilpotent variables. This implies that $(-B) as a function of these variables is a finite 
power series. However, when replacing B with p we continue this power series to an 
infinite one. It is important to realize that this step is not problematic at all, because 
$(5) = ^{K) results from the Fourier transform of P{H). To illustrate the feasibility 
of this continuation, we refer to the cases covered by Eq. ()H1|1 . 
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When writing out the expression trgcrp, one sees that the imaginary units due to 
the Wick-type-of rotation in the fermion-fermion blocks nicely combine to —1 such 
that the whole expression trg ap is real. This is of course necessary to make the Fourier 
transform well-defined. The inverse of the Fourier transform ()36|) reads 



We notice that the prefactor 2^^*^'^"^^ does not involve vr, because we have the same 
number of commuting and anticommuting variables. Due to the invariance of the 
measure d[p], the rotation invariance of the characteristic function $(p) gives with 
Eq. ()37|) directly the same property for the probability density, such that 



The rotation invariance of P{H) implies the corresponding feature for Q{cr). 

There is a good reason why we defined Q(o") as above. Nevertheless, what we need 
now to carry through our construction, is the integral representation 



of the characteristic function considered as a function of B. At first sight, there is a 
problem, because trg aB is not real anymore. As the imaginary unit is present in the 
fermion-fermion block of a, but absent in that of B, the Fourier integral (jH^ seems ill- 
defined. However, as argued above, one can also Wick-rotate the relevant elements of 
B. Even if one chooses not to do that, everything is under control, because the matrix 
elements CpCg of B in the fermion-fermion block are nilpotent. The corresponding 
expressions in exp (ztrg aB) consist of a finite number of terms, and no convergence 
problem for the a integration can occur. 

3.5. Generalized Transformation Formula 

After these preparations, we are in the position to derive a transformation formula 
which expresses the probability density in superspace as an integral over the probability 
density in ordinary space. Using the result pHj) . we have 




(37) 




(38) 




(39) 



N 



tl HY= Y^HnnYn 



n=l 
k 



k 



p=i 



{iH(k+p)ik+p)) {iyp2) = trghy 
p=i 



(40) 



with 



h = diag {H] 
Hence, Eq. fITBj) yields 



, Hkk, ^-f^(A:+l)(fc+l), • • • , iH{2k){2k)) ■ 



(41) 



$(5) = ^K) 



c?[i7] P (if) exp (ztrg hy) 



(42) 
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A proper definition of tlie diagonal matrix h made it possible to employ the supertrace 
in Eqs. fl40|) and ()42|) . We plug the characteristic function into the inverse Fourier 
transform ()37|1 and find 

Q{a) = 22^(^-1) J d[p] exp (-itrg ap) J d[H] P{H) exp (itrg hr) , (43) 

where we use that y and r have precisely the same form. Assuming that the order of 
integrations may be interchanged, we arrive at the generalized transformation formula 

Qia)= j d[H]P{H)x{cy,h) . (44) 

The function 

X((T, h) = 22^(^-1) J d[p] exp (itrg {hr - ap)) . (45) 

is a projector which is related to, but different from a 6 function. It might look surprising 
that the integrand contains the full matrix p as well as its eigenvalue matrix r, but 
recalling the derivation, this is rather natural. The term exp(itrg/ir) stems from the 
Fourier transform of the probability density P{H) in ordinary space. Although it is 
conveniently written in a supersymmetric notation, it is exclusively rooted in ordinary 
space. Thus, anticommuting variables may only implicitly be present, which makes it 
plausible that r appears, but not the full p. 

The projector satisfies the important normalization property 

J d[a] x(a, h) = J d[p] 5^'^'^ (p) exp (ztrg hr) = I , (46) 

where 5^'^^^\p) is the product of the 5 functions of all Ak"^ independent matrix elements 
in the supermatrix p. This then gives directly the normalization 

j d[a] Q{a) = J d[H] P{H) J d[a] x(a, h) = J d[H] P{H) = 1 (47) 

of the probability density in superspace. As one should expect, the normalization of 
P{H) yields the normalization of Q{a). One is tempted to conclude that this feature 
wraps up the whole convergence discussion if the choice (fH^ has been made. Such an 
interpretation is corroborated by the character of the projector x((T, h) which is truly 
convergence friendly under the integral. Nevertheless, this thinking comes to terms 
when considering the complexity of all intermediate steps. Unfortunately, it prevents 
us at present from providing the impression stated above with more mathematical 
substance for a general P{H). One possible problem is related to the Wick-type-of 
rotation. All invariants ii = tr E'^"^ are positive semi-definite for all integer m. 
This is clearly not so for the corresponding invariants trgcr^'" = trgs^™', where we have 
positive semi-definiteness only for odd integers m. This does not inevitably lead to 
difficulties, because a term exp(— trif^m) P{H) is not necessarily mapped onto its 
analog exp(— trg cr^™) in Q{a), but it illustrates at which points problems could arise. 
Nevertheless, anticipating the discussion to follow in Sections 0] and |S1 we mention 
already now that the whole problem can be considered exclusively in Fourier superspace 
such that only the convergence properties of the characteristic function matter. 
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3.6. Generating Function 

Having obtained the probability density we use Eqs. (l30|) and p9|) in formula (|TT|l . 

The remaining steps to be done are then exactly as in Ref. and we arrive at the 
result 

Zkix + J) = j d[a] g(a)detg (a - - j) , (48) 

where 

Q{cr)= j d[H]P{H)x{cr,h) (49) 
is the probability density in superspace. 

3. 7. Norm- dependent Ensembles Revisited 



The transformation formula (j49j) generalizes a transformation formula which we obtained 
for norm-dependent random matrix ensembles PU] • We revisit this case to acquire some 
experience with the generalized transformation formula. The probability density P{H) 
of a norm-dependent ensemble depends on H only via tvH^. In Ref. fSl, the class of 
these ensembles was constructed by averaging Gaussian probability densities over the 
variance t, 

oo 

PiH) = J f{t)^^^j^-^,e^V (-^tr/f^) dt . (50) 

where the choice of the spread function f{t) determines the ensemble. With the 
transformation formula (jlHI), we find 

oo 

g(a) =22^-('=-i) j dtfit) J d[p] exp{-ttTgap) 



/ d[H]^^j^^^^exp (^-Itrif^j exp(ztrg/ir) 

oo 

^22kik-i) J ^^j^^^ J ^[^j exp(-itrgap) 



oo 

= 2^(^-^) J dtfit) j d[p\ exp(-itrgap) 2^('=-i) exp (^-^trgr^^ 



oo 

= |rft/(t)2'=('=-i)exp(--trga2) (51) 



which is indeed the correct result. We mention in passing that it allows one to express 
the mapping of norm-dependent ensembles from ordinary to superspace as one single 
integral in terms of the probability density alone [Tpj, i.e. without explicit appearance 
of the spread function. 
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4. Supersymmetric Formulation in Fourier Superspace 

Another supersymmetric formulation of the generating function will prove most helpful 
for calculations of the correlation functions later on. Also from a conceptual viewpoint, 
it has some rather appealing features. In Section 14. H we construct the new formulation 
by exploiting a convolution theorem, and in Section 14.21 we give a direct derivation. 



4-1- Applying a Convolution Theorem 

According to Eq. (jlH|l . Zk{x + J) is a convolution in supermatrix space. For three 2k x 2k 
Hermitean supermatrices a, p, t and for two well-behaved functions Qiic), g2{<y) as well 
as their Fourier transforms Gi{p), 6*2 (p), one easily derives the convolution theorem 

I d[a] gi{a)g2{r - a) = 2^'^'-'^ J d[p] exp (-^trg rp) Gi(p)G2(p) . (52) 

In the present case, we have t = x + J. We already know the Fourier transform of Q{o'), 
it is just the characteristic function $(p). The Fourier transform 

/(p) = J d[a] exp (itrg per) detg ^^a^ . (53) 

of the superdeterminant is needed. It can be viewed as a supersymmetric generalization 
of the Ingham-Siegel integral, whose ordinary version has recently been used in the 
framework of supersymmetric methods JE]. Obviously, /(p) only depends on the 
eigenvalues r of p. In [Appendix C[ we show that it is given by 

/(p) = CNk n Q{rpi){irpi)^ exp (-er^i) ^ i^^f"^ 

p=l ^^p2 

1 /z27r(-l)^-i\' , , 

CNk - ^^(^ 1^ (iV-1)! ) ■ ^^^^ 

We notice that /(p) is almost equal to detg'^'^p, apart from the restriction to negative 
eigenvalues rpi and the occurrence of the functions 5(rp2) instead of 1/r^. Loosely 
speaking, the Fourier transform maps the superdeterminant raised to the power —N 
onto the superdeterminant raised to the power +A^. We find from Eqs. (j48j) and (j52j) 

Zk{x + J) = 22'=(^-i) I d[p] exp (-ztrg {x + J)p) ^{p)I{p) . (55) 

Thus we arrive at the remarkable insight that only the characteristic function (^{K) 
is needed in the generating function and, thus, for the calculation of the correlation 
functions. It is of considerable conceptual interest that the probability density in 
superspace Q(cr) follows in a unique way from the one in ordinary space P{H), but the 
use of Q(cr) can be avoided if the Fourier superspace representation is more convenient 
in a particular application. 
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4-2. Direct Derivation 

Since Q{cr) does not appear anymore in the expression (j55|l . the question arises if one 
can obtain this result without going through the construction of the probabihty density. 
This is indeed possible. We go back to Eq. (^TJ) and do the average over the ensemble, 

Zk{x + J) = Yl / d[zp] exp (iLpzliiLpE - Xp + Jp)zj)j 
p=i'' 

J d[Cp] exp {tClizLpB -Xp- Jp)Cp) ^{K) . (56) 
We now use the insights of Section 13.41 and insert an integral over a 5 function, 
^{K) = $(5) 

= Jd[p]^p)6^'''\p-B) 

= 22^('=-i) J d[p] ^p) J d[a] exp (-ztrg a{p - B)) , (57) 

where p and a are 2k x 2k Hermitean supermatrices, to which the Wick-type-of rotation 
has been applied in the fermion-fermion blocks. Again, one might argue that this makes 
the expressions in Eq. ()57p ill-defined, because these matrices and the matrix B are 
treated on equal footing, although no Wick-type-of rotation has been applied to the 
latter. The same reasoning as in Section 13.41 can be employed: Either one also Wick- 
rotates B or one argues that the integrals in Eq. (|37j) are well-defined because the 
elements of B in the fermion-fermion block are in any case nilpotent. We plug Eq. (jFTj) 
into Eq. The integrals over the vectors Zp and (p can then be done in the usual 

way, and we have 

Zk{x + J) = 22'=('=-i) J d[p] $(p) 

J d[a] exp (— itrgap) detg~^ (^a — x^ — 
= 22'=(*^-i) J d[p] $(p) exp (-itrg (x + J)p) 

J d[a] exp (— itrg ap) detg^^a^ , (58) 

where we shifted a by a; + J in the last step. The remaining a integral is, after changing 
variables from a to —a, precisely of the Ingham-Siegel type and we obtain Eq. 

Of course, the probability density Q(o") is somewhat hidden in Eq. (jHTj) . However, 
to actually obtain it, one has to do the p integral, which would require an interchange 
with the a integration. Avoiding the introduction of the probability density Q(cr) in the 
derivation sheds new light on the convergence issues. If P{H) is a Schwartz function, 
^{K) is a Schwartz function as well and the convergence discussion can be exclusively 
restricted to the Fourier superspace and to the properties of the characteristic function 
when passing from ordinary space, i.e. from ^{K), to superspace, i.e. to $(-B) and $(p). 
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5. Correlation Functions in Terms of Eigenvalue Integrals 

In Section 15.11 we briefly review the integrals that we need over the unitary group in 
ordinary and in superspace. We derive a first general result by identifying fundamental 
correlations in Section 15.21 In Section 15. 3| we carry out the procedure outlined 
in Section 12.31 and obtain supersymmetric integral representations of the correlation 
functions for arbitrary positions of the imaginary increments. Another general result 
is given in Section 15. 4^ exclusively in terms of eigenvalue integrals. In Section 15. 5^ we 
discuss a probability density involving higher order traces as an example. 



5.1. Eigenvalue-angle Coordinates and Group Integrals 



The Hermitean random matrix is diagonalized according to H = UEW with E = 
diag {El, . . . , En) and with U being in SU(A^). The volume element in these coordinates 
reads 

iV(7V-l)/2 

d[H] = A%{E)d[E]dfi{U) , (59) 

' • 1 ln=l 

where we introduced the Vandermonde determinant 



Ajv(^) = det 



E. 



m—l 



n,m=l,...,N 



n {En - E, 



(60) 



n<r?i 



The invariant measure dii{U) in Eq. (|59|) is normalized to unity. We will make use of 
the Harish-Chandra-Itzykson-Zuber integral |IT71 



/ rrprrtp^ Hl'/n! det [exp(^E„i?^)]„ ^ 



(61) 



■ ,m=l,...,Af 

where R = diag {Ri, . . . , R^) is also a diagonal matrix. In particular, we will employ 
this result for the case that i?„ = for n > 2k. This can be obtained in various ways, 
as for example in Ref. |19j . 

N-l 

n 

n=N-2k+l 

exp{tEr.R2k) 1 ■ ■ ■ E^-^'^-^ 



dn{U) exp{itTUEU^R 



det 



exp{iEnRi 



n=l,...,7V 



(62) 



AN{E)A,,{R)Uli,R^-'' 

where we write R = diag (i?i, . . . , R2k)- 

In superspace, the diagonalizations of the Hermitean supermatrices a = usu^ and 
p = vrv^ have already been introduced in Eqs. (jH^ and (jMj) . The volume element d[p] 
reads in eigenvalue-angle coordinates ^2] 

d[p] = Bl{r)d[r]dp{v) , (63) 

where the function 



BJr) 



Afc(ri)Afc(zr2) 
Ilp<q{rpi - ir g2) 



det 



1 



(64) 
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is the superspace equivalent of the Vandermonde determinant. The supersymmetric 
analog fT^ 120] of the Harish-Chandra-Itzykson-Zuber integral is given by 

dn{v) exp (jtig vrv'^s 



det fexnf?'r^o,s„oll 

(65) 



det [exp(zrpiSgi)]p_^^^^ ^^det [exp(irp2S52)]p,g=i 



As in Refs. IT^ . the normalization of the invariant measure dfi{v) is chosen such 
that formula when applied to a shifted Gaussian distribution, yields the proper S 
function in the curved space of the eigenvalues for vanishing variance. 

5.2. General Result as an Average over the Fundamental Correlations 

The supergroup integral can now directly be applied to the Fourier superspace 
formulation (j^ . because both of the functions $(p) and I{p) depend only on the 
eigenvalues r. This is the merit compared to the original superspace formulation PHj). to 
which the result cannot be applied in general. In the case of a Gaussian probability 
density, a shift of the integration matrix o" by a; + J gives a form suited for the application 
of the supergroup integral [121 CSl El ■ In the general case, however, Eq. is much 
more convenient. We find 



Zk{x + J) = 1 



2fc(fc-i) / ^ N 



Bk{x + J) \2tt, 

J (i[r] i?fc(r) exp (—itrg (x + J)r) $(r)/(r) . (66) 

Two remarks are in order. The first term, i.e. unity, stems from a certain boundary 
contribution which only appears in superspace. In physics, it is often referred to as 
Efetov-Wegner-Parisi-Sourlas term ^1 1^ |221 12^] , while it goes by the name Rothstein 
contribution |21] in mathematics. In the present case, it yields the normalization 
Zkix) = 1 of the generating function, because one easily sees that 1/ Bi^i^x + J) vanishes 
at J = 0. Formally, the boundary contribution is obtained by putting p = in the 
integral (j^ . by using $(0) = 1 and /(O) = l/2'^(*^~^) according to Eq. and 
to Appendix C[ and by finally dividing the result with the factor 2^^^~^^ which is due to 



our definition of the volume element d[p]. There are various methods to explicitly justify 
this procedure in the case k = 1. In Ref. [221, for example, it is directly constructed 
from Rothstein's theorem. However, there is a problem, because none of those explicit 
methods could be extended so far to our eigenvalue-angle coordinates for k > 1. We 
can thus not exclude that further boundary contributions exist. Nevertheless, as to be 
discussed below, we are confident that they are not important for our purposes. 

The second remark concerns the determinants in the formula ()65|) which are not 
present in Eq. As the functions $(r) and /(r) are invariant under permutations of 

the variables rpi as well as of the variables rp2, it suffices to keep only one term of each 
determinant, because all others yield the same under the integral. 
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We can now proceed in different ways. Here, we begin with inserting the 
characteristic function in the form 

$(r) = J d[H] P{H) exp (ztrg hr) (67) 

as given in Section I3.51 Upon interchanging the r and the H integral we find the 
expression 

Z,{x + J) = 1 + J'^^l^j, I d[H] P{H)Ut''\x +J-h), (68) 

where we introduced the fundamental correlation function 

'^'i)(s) = 2'=('=-i) J d[r] Bk{r) exp {-ttig sr) I{r) (69) 

as a new object. In Eq. (|HHjl . we have to set s = x + J — h. We refer to the 
correlation function (jU^ as fundamental, for it gives all structural information about 
the correlations before averaging over the probability density P{H). The fundamental 
correlation function is the Fourier transform of the function /(r) in the curved eigenvalue 
space. It is closely related to the backtransform of /(p), i.e. to the superdeterminant 
detg'^cr^, but it is not quite the same. We discuss that in [Appendix D[ 

The result ()68|) is not a trivial reformulation of Eq. defining the generating 
function. While it is obvious from Eq. that only the eigenvalues of H are relevant 
for the ensemble average, Eq. (jUHj) makes a different statement, namely that only 2k 
diagonal elements of H enter the computation of the average. Using the determinant 
structure ()64j) of Bk{r) and formula ()54|1 . we find that the fundamental correlation 
function has the determinant structure 



i?™(s) = det 



J (fund) 

where the fundamental kernel is given by 



C^'--''\s,,,ts,2)] ^ , , (70) 



e(n)(-i)"^i^^ . (71) 



We suppress the indices p and q in the integration variables ri and r2- It is shown 
Appendix D that the fundamental kernel can be written as 



m 



1 N—l /; yi „ 

C^'''^''\s,,,is,2) = - - E / dr,{ir,r e^v i-^rlS-l 



^ ^0 i-^pi) 



As this is a finite geometric series, we may also write 

^(fund). ^ ^S2)= ^ ^~(^^^2/v)^ ^ 1 (a^i)^ - {tSg2)^ 

TTSpi l-{iSg2/Spi) 7r(spi)^ s-^-iSg2 
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The fact that the fundamental kernel has a representation as a finite series and as a 
ratio of differences is reminiscent of and related to the Christ off el-Darboux formula [22] 
in the theory of orthogonal polynomials. 

The correlation functions according to Eq. Q are then quickly obtained using the 
steps of Ref. [12]. We find 

Rk{x,, ...,Xk)= I d[H] PiH)Rt'"'\x - h) 
' d[h] P^'''''\h)Rf''"'\x - h) 



det[C^'''^''\xp-Hpp,x,-iH^k+,)i^u+,))\^^^^ ^. (74) 

The correlation functions are convolutions of the the fundamental correlations with the 
reduced probability density 

pi^--^)[h) = j d[H/h]P{H) (75) 

found by integrating P{H) over all variables except the 2k diagonal elements h of H. 
The result ()74|1 holds for arbitrary rotation invariant Random Matrix Ensembles. We 
notice that the reduced probability density is connected to the characteristic function. 
One sees that either directly from Eq. (jFTTj) or by performing the following steps, 

p(rcd)(^^ = I d[H'] 6{h' - h)P{H') 

= J d[r] J d[H'] exp (^trg {h' - h)r) P{H') 

= J ^^'^^ {-itighr) <l>(r) . (76) 

Hence, P^^'^^\h) is the Fourier backtransform of the characteristic function depending 
on the 2k coordinates r — which are here viewed as describing a fiat space — onto 
a function defined in the fiat space with coordinates h. This is very different from 
Fourier transforms in curved space, when the eigenvalues r are interpreted as the radial 
coordinates of a Hermitean supermatrix p. 

It is somewhat surprising that the probability densities in ordinary space P{H) 
or P^^'^'^\h), respectively, suffice to write down Eq. One might conclude that 

this obliterates the above convergence discussion related to the functional forms of the 
probability densities in ordinary and superspace and of the characteristic function. 
Formula ()74|) indeed gives reason to be optimistic. However, we recall that the 
characteristic function was used in the derivation, even though it does not appear any 
more explicitly. Hence, we must require the existence of the characteristic function and 
also that the H and the r integrations can be interchanged when going from Eq. ([67jl 
to Eq. 

The inherent determinant structure ()70p of the fundamental correlations will be 
destroyed in general when averaging over the Random Matrix Ensemble. It will be 
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preserved if the reduced probability density factorizes according to 

2k 

p^'^^^^h) = n p^'^'^^H^pp) • (77) 

p=l 

However, this is not the only situation in which the determinant structure survives. The 
Mehta-Mahoux theorem implies that the correlation functions Rk{xi, . . . ,Xk) can be 
written as determinants for all rotation invariant probability densities which factorize 
in their eigenvalue dependence, 

TV 

PiH) = P{E) = II P^'^^^En) . (78) 

n=l 

One would not expect that the the factorizations (f77|) or (f7S|) are completely 
independent, but we have not looked into this further. In the present context, it is more 
important that the applicability of the Mehta-Mahoux theorem is limited to precisely 
the case when the factorization (f7H|) holds. It is thus a quite attractive feature of the 
result (f7^ that it is valid for all rotation invariant probability densities which have the 
property P{H) = P{E), but which do not need to have any factorization property as in 
Eqs. fl77|) or (fTSjl . In this sense, formula (|7^ is more general than the Mehta-Mahoux 
theorem. In Section E31 we will give another result, also valid for all rotation invariant 
probability densities. Since it is formulated in terms of integrals over the eigenvalues, 
its structure is somewhat different from formula ()74j). 

As an easy check of our findings, we show in [Appendix E that Eq. ()74j) yields 
immediately the GUE correlation functions. This is important, because it strengthens 
our confidence that we treated the Efetov-Wegner-Parisi-Sourlas term EH 122 ISH] 
consistently. 



5.3. Correlations Functions of the Imaginary Parts and for Arbitrary Positions of the 
Imaginary Increments 

As discussed in Section ITTl the correlation functions /^^(xi, . . . , x^) as defined in Eq. (jT)) 
are the main object of our interest. We now construct integral representations for them 
and, in addition, also for all correlation functions Rk(xi, . . . ,Xk) as defined in Eq. © 
for arbitrary positions of the imaginary increments. To avoid introduction of hyperbolic 
symmetry, we restricted ourselves from Section IH.4I on to the case that all imaginary 
increments lie on the same side of the real axis. However, applying the strategy outlined 
in Section |2ini we can recover every correlation function that we want. 

It is convenient to use the general result (fTIj). allowing us to conduct the 
construction by only looking at the fundamental correlation function Rf^^^\x — h). 
Due to its determinant structure, it depends on one fixed energy Xp either in 
the form C^^^^'^^Xp — Hpp,Xp — iH(^k+p){k+p)) or in the form C'^^™'^)(xg — Hgg,Xp — 
^i^(.+p)(.+p))^(^^"'^Hxp - Hpp, X,, - where g ^ p and g ^ V- From^the first 

of the expressions ([7^ we conclude that in both cases the dependence of r!>^^^'^^ [x — h) 
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on the fixed energy Xp is a finite sum of the terms 

oo 

^nmixp) = (iH(^k+p){k+p) - Xp^ J (iri(iri)"exp (-iri{Xp - Hpp)^ , (79) 



where n = m is possible. As the average over the ensemble is linear, it suffices 
to investigate the functions Anm{xp) in order to study the energy dependence of the 
correlation functions Rk{xi, . . . ,Xk). According to Section 1221 we study the Fourier 
transform 

^ +00 

^nm(tp) ~ 7^= / dXp GXp (^itpXp) Afimi^p) ■ (^0) 

\' In J 

— oo 

Shifting Xp by Hpp, it can be cast into the form 

^ +00 

p) — —7= exp (itpHpp) / dxp exp (itpXp) 

\/ 2n J 

—00 

00 

{iH(k+p){k+p) - Hpp - Xp) j driiiriY exp (-irix^ ) 



1 / (9 \ ™ 

= 6xp (itpHpp) iH(^i,+p){k+p) — Hpp + i— 
2tc \ dtp j 

CO +00 



(iri(iri)"' exp (— eri) / c?Xp exp (ixp(tp — ri)) 



m 



27rexp [itpHpp) ^iHi^k+p){k+p) - Hpp + i-^ 



dri{iri)'^ exp {—eri) 6(tp — ri) . 



As the ri integration extends over the positive real axis only, the integral is zero whenever 
tp < 0. All derivatives are zero as well in this case, implying that the entire expression 
is proportional to Q{tp) . For tp > 0, the integral yields (itp)" exp(— etp). All derivatives 
of the exponential function give terms containing powers of e and thus vanish in the 
limit £ — >■ 0. Here, we may assume that the tp integral cannot yield bare singularities 
in e. We can thus neglect all these terms and write exp{—etp) in front of the entire 
expression. We find 

Anm(tp) = i2Q(tp) exp {-etp) \nm{tp) (82) 

where 

"\/27r ( (9 \ ™ 

^nm{tp) = — ^exp [itpHpp) y,H(_k+p){k+p) — Hpp + ig^j {^^pY ■ (83) 

Indeed, Eq. (jH^ directly implies expression ^ and we can read off the desired integral 
representations. The function Xnmiip) is recognized as Fourier transform of 

+ 00 

Km{.Xp) = - i {iH(^k+p)(^k+p) - Xp)"" J driiiri)"" exp {-iri{xp - Hpp)) , 
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Qn 



-l)"27r {iH^k+p){k+p) - Xp) -jr^Sixp - Hpp) 



, m _ 



= - 2n\ {iH^k+p)ik+p) - Xp)" Im _ . (84) 

Collecting everything, we arrive at 

...,Xk)=J d[h] P^'''''\h)Rt'"'\x - h) (85) 
with the fundamental correlation function 



i?™(s) = det 



-.(fund) 

and the fundamental kernel 



C(f"°<^)(spi,^^.2)] ^ , (86) 



1 """^ 



C^''''"'\sp,,is,2) = -T^E^^ / rfri(^ri)"exp (-zri.pl) 



... 1 



-EK2)"Im— — ^. (87) 

n=0 v-^pl/ 

Hence one simply has to replace the singularities l/(s~i)"^^ everywhere with their 
imaginary parts. Tracing back these considerations, we realize that all necessary 
modifications reside in the r^i integrals and specifically in the function /(r). Replacing 
Eq. (jni) with 

is equivalent to the above discussed steps made to obtain Rk{xi, . . . , Xk). 

With the help of formula ()lUp. it is now an easy exercise to construct integral 
representations for the correlation functions Rk{xi, . . . ,Xk) defined in Eq. (j2)) with 
arbitrary positions of the imaginary increments. Formulae (ffn|) and (ff^ remain valid if 
the fundamental kernel is replaced with 

1 {^s,,T °° 



^('"'^'^^^,^^.2) = T - E / ^ri(2ri)'^exp {T^r^s 

n=0 ■ n 



where the notation s^i indicates that the imaginary increment is chosen according to 
x^ = Xp — iLpE. In terms of the function /(r), this is equivalent to replacing Eq. (j54j) 
with 

Hp) = CNk n Q{Lprpi){irpi)^ exp {-LpErpi) — . (90) 

p=l '^^p2 

Thus, we obtain supersymmetric integral representations for all these correlation 
functions without using hyperbolic symmetry. 

The insights just presented may be viewed as a more formal justification of the 
procedure denoted by the operator symbol ^5 in Refs. [121 HBl IHj- We argued in 



Arbitrary Rotation Invariant Random Matrix Ensembles and Supersymmetry 



25 



these studies that the generating functions satisfy a diffusion process. The diffusion 
propagator contains no information about the positions of the imaginary increments, 
this is exclusively contained in the initial condition of the diffusion. Moreover, the 
diffusion propagator is nothing but the supersymmetric Harish-Chandra-Itzykson- 
Zuber integral over the unitary supergroup, not involving any non-compact degrees 
of freedom. This can be verified in an elementary way by simply plugging it into 
the diffusion equation. Hence, one is free to adjust the positions of the imaginary 
increments as needed, which essentially defined the operator 53. We have now given 
another justification. Nevertheless, it remains an interesting mathematical question to 
also derive all that from group integrals involving non-compact degrees of freedom. 



5.4- General Result in Terms of Eigenvalue Integrals 



A further integral representation follows from Eq. (jbbj) . We take the derivatives with 
respect to the source variables as in Ref. ^2] and in Section and find 

Rk{xi, . . . , Xk) =2^^^-^^ f d[r]Bkir)exp{-itigxr)^{r)I{r) . (91) 



The correlation functions Rk{xi, . . . ,Xk) as well as those for arbitrary positions of the 
imaginary increments are obtained as in the previous Section 15.31 we simply have to 
replace J(r) according to Eq. (jK^ by /(r) according to Eqs. (jHHjl or (jnUjl . respectively. 
We expand the determinant Bk{r) by introducing the permutations u of the indices 
p = 1, . . . , k and write 

i?fc(xi,...,a;,)=2'=('=-^WE(-l)'^"^ 

d[r]^{r)[[ ^ ^ 



rpi - tr^^p)2 



Q{Lprp^){irp^)'' exp (-L,er,i) '^^^^^ , (92) 

where j{uj) is the parity of the permutation uj. The 5 functions allow us to do the k 
integrals over the variables immediately. We integrate by parts and use Leibnitz' 
rule to work out the derivatives of products, 

i?fc(xi,...,a;fc) = («27r)'=E(-l)^('^) 



1 r 

— -I drii e(Lirii)exp (-zxirii - LiCTii) (-irii)"i 

N-l ^ +9° 

V — - I drkiQ{Lkrki)exp{-ixkrki - Lkerki) 

(93) 



-^rfci)"'= 



-oo 



dK^i exp (- Ep=i a:^(p)rp2) $(r) 



r2=0 
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This result is valid for an arbitrary rotation invariant probability density. The structure 
of this expression is quite different from the one in Section 15.21 where the correlation 
functions were found to be a convolution of the reduced probability density with the 
fundamental correlations. 

It is instructive to see how the correlation functions can acquire a determinant 
structure. An obvious feature leading to this would be a factorization 

k 

$(r) = n$(-)(rpi)<l>(^^)(rp2) (94) 
p=i 

of the characteristic function. We find immediately 



Rk{xi, . . . ,Xk) = det C{xp,Xq) (95) 

L Jp,5=l,...,fc 



with the kernel 



Af-l 



C(^Xp, Xq) y ] 



1 9"exp(-Xgr2)$(''^)(^2) 



vr n\ drl 



r2=0 



J dri e(Lpri) exp {-ixpr^ - LpSri) (-zri)"$(^")(ri) , (96) 

— oo 

where we suppress the indices p and q in the r variables. We notice that the GUE 
case is trivially recovered. We then have $^''^-*(rpj) = exp(— rp^/4) which combines 
in the derivative expression with the exponential to the generating function of the 
Hermite polynomials, and the integral yields the generalized Hermite functions as 
given in [Appendix E| It is conceivable that mechanisms other than following from the 
factorization (|94p can be identified that also lead to a determinant structure. However, 
as the merit of Eq. is its completely general character and its independence of such 
factorizations and determinant structures, we have not explored this issue further. 

One can wonder whether it is helpful to integrate over the group SU(A^), i.e. over the 
ordinary unitary matrix U diagonalizing H, before inserting the characteristic function 
$(r) in formula With the help of Eq. ^ we find 

^N{N-l)/2 



Hr) = ^, ^^„i ^, J d[E]Al{E)P{E) J df,{U)exp {HtUEU^r) 

^NiN-l)/2 



iiN-mk-i)N\ Yi^-^^ n\ A2,(ri, r^) nf=i(vr,2)^-2'= 
d[E]AN{E)P{E) 
exp(zE„rn) ■■■ exp(zKrfc2) I ■ ■ ■ E^-^^-"-] _^ (97) 



det 



where we have to set Rp = rpi, Rp+k = rp2, p = 1, . . . ,k. The eigenvalues rp2 do not 
come with an imaginary increment in the formula above. This is also important in 
^2k{ri,r2) which is the ordinary Vandermonde determinant of the 2k variables ri and 
r2. As the whole integrand is invariant under permutations of the eigenvalues En, we 
may replace the determinant stemming from the group integration by the product of 



Arbitrary Rotation Invariant Random Matrix Ensembles and Supersymmetry 



27 



its diagonal elements, because all other terms yield the same result. However, from the 
resulting expression 

^N{N~l)/2 



^(r) 



d[E]A^{E)P{E) 

k N 

l[exp{tiE,r,, + Ek+prp2)) U E^-^^-^ (98) 

p=l n=2fc+l 

it is not immediately obvious anymore that its limit of vanishing rpj remains finite, given 
by the normalization $(0) = 1. 

We give the expressions ()97|) and ()98p mainly for the sake of completeness, because 
they are not particularly useful in their general form. Although the powers in the 
denominator are not real singularities in Eqs. ()97p and (j^Hjl . they become truly singular, 
if one tries to exchange the order of integrations and to do the r integrations first. 

5.5. Ensembles Involving Higher Order Traces as an Example 

As it might be helpful to illustrate our findings by an example, we consider the 
probability density 

P{H) = bM,M, {tTH^')''\xp (-tiH^) (99) 

for a fixed pair of integers Mi,M2 = 0,1,2,.... The constant &AfiAf2 ensures 
normalization. The Gaussian case is recovered for Mi = or M2 = 0. A few obviously 
meaningless cases have to be excluded, such as the choice Mi = M2 = 1, which makes 
the normalization integral vanish. While the probability density (j^n]) is still in the 
norm-dependent class discussed in Ref. ^U] for Mi = 2, it is not for other values of 
Ml. Importantly, the probabihty density (j99|l does not factorize according to Eqs. (ffTjl 
or (fTHj) . In particular, this means that this Random Matrix Ensemble is not covered 
by the Mehta-Mahoux theorem, although we do not exclude that is possible with some 
efforts to extend the latter properly. Formula ()74|1 provides a direct way to calculate 
the correlation functions for such ensembles. However, as we aim at addressing the 
conceptual issues in the present contribution, we refrain from presenting the quite 
cumbersome expressions too explicitly. We rather sketch the calculation briefly and 
infer what kind of structure the correlations functions will acquire. It is obvious from 
definition (ffH| that the reduced probability density has the form 

p('-^'i)(/i) = exp (-trg h^) J2 a{m}S{m}ih) , (100) 

{m} 

where the a^m} are constants and where 

2k 

V}W=En^p^"''' (101) 

UJ p=l 

are symmetric functions, i.e. linear combinations of products involving a set {m} of 
integer exponents rup, symmetrized by summing over all permutations lj of the indices 
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p = 1, . . . , 2k. The exponents rup are between zero and Mi + M2 with the restriction 
that their sum does not exceed Mi + M2. Thus, the correlation functions are given by 

Rk{xi, . . . , Xk) = XI '^W / ^i^l exp ti h"^^ S{m}{h)Rk'''"^\x-h) .(102) 

{m} 

Upon inserting Eq. pOlj) and using the determinant structure of the fundamental 
correlations, we obtain 



^fc(xi,...,a;fc) = ^ai^l^det p^^(^)^^(,^^)(xp,a;g)l (103) 

{m} ^ Jp,<Z-i,...,fc 



where the kernel 

1 1 

Cmim2 ("^p? -^g) exp ^ "''p) ^ ^ j'^nmi (•^p)^nm2 ("''g) (-^0^) 



has a structure formally similar to that of the GUE kernel. The functions 

+00 00 

Vnmiixp) = J dHppexjp (^-Hpp^ H^^ J dri{ir i)"" exp (^Tir lixj - Hpp 

-00 
+00 



rm2 

k+q)(k+g)J ^^{k+q)ik+q) 



^Xq - iH(^k+q){k+g)) (105) 

can be written as finite weighted sums of the generalized Hermite functions which 
are discussed in [Appendix E and of the ordinary Hermite polynomials, respectively. 
According to the result ()103p . the correlation functions are linear combinations of 
determinants. 

Alternatively, this calculation can be carried out using the results of Section It 
follows from the inverse of formula (fTBjl that the characteristic function has a form very 
similar to the reduced probability density, 

$(r) = exp (-^trgr^] ^ (106) 
^ ^ ^ {m} 

with new constants ajm}- With the help of Eq. or (jH!^ this leads in a straightforward 
manner to the above mentioned linear combinations of determinants. 



6. Summary and Conclusions 

We derived supersymmetric formulations for arbitrary rotation invariant Random 
Matrix Ensembles. The construction is based on an algebraic duality between 
ordinary and superspace which made it possible to generalize the Hubbard-Stratonovich 
transformation. We identified an integral transformation that involves a projector and 
yields the probability density in superspace from the one in ordinary space. However, 
we showed that despite the conceptual insights thereby obtained, the theory can be 
formulated without using the probability density in superspace. It turned out that it 
is possible and often even better to work in Fourier space, because the characteristic 
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functions of the probability density have the same functional form in ordinary and 
superspace. At present, it appears to us that it is a priori easier to analyze some 
convergence issues in Fourier superspace, but to make more definite statements will 
require additional work. It is not inconceivable, that manifestly invariant theories can 
be constructed in Fourier superspace. 

This leads us to the symmetry issue. There seems to be no way around hyperbolic 
symmetry if one wishes to set up non-linear a models. Here, however, we were interested 
in exact, non-asymptotic results. Although this requires the calculation of certain group 
integrals, it simplifies the symmetry: We showed that compact supergroups are sufficient 
to construct supersymmetric integral representations of the correlation functions for 
arbitrary positions of the imaginary increments. This is a more formal justification of 
a procedure which we have been using in previous work. We conclude that hyperbolic 
symmetry is a necessity for non-linear a models, but not for supersymmetric theories in 
general. Nevertheless, even though mathematics can be nicer than one expects, it is an 
interesting challenge to also derive those supersymmetric integral representations from 
a version of the theory in terms of non-compact groups. 

We gave two general results for the correlation functions. The first one involves 
certain correlations to which we refer as fundamental, while the second one is only in 
terms of eigenvalue integrals. Both results are valid for arbitrary rotation invariant 
Random Matrix Ensembles. In particular, no factorization property of the probability 
density has to be assumed. 
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Appendix A. Equality of the Traces 

The assertion (jUj) is obviously correct for m = 1, because we have 

k k 

tiK = ^ tr (^LpZpzl - CpCj) = {Lpzlzp + CjlCp) 
p=i p=i 

= tTgB . (A.l) 

For m = 2, 3, . . ., we find 

tr K"" = tr ALA^ ■ ■ ■ ALA^ 

= tr AL'/^L'/^A^AL'/^ ■ ■ ■ L'/^A^AL'/^L'^'A^ 

= iiAL^'^B'^-^L^'^A^ . (A.2) 

Without anticommuting variables, we could now simply use the invariance of the trace 
under cyclic permutation and would arrive at the desired result (j221), but with an 
ordinary trace also on the right hand side. To carefully account for the anticommuting 
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variables, we write C = 5'"^^ and introduce the upper indices (cl) and (c2) in boson- 
fermion block notation for the commuting variables as well as (al2) and (a21) for the 
anticommuting ones. We obtain 

tri^- = Etr {z,Ll/^CifLl/^zl - z^Ll/^C^f^Cl 

, /• ^(a21) rl/2 t _ ^ (-<{c2)^f\ 
' Sp^pq ^q Sp^ pq Sq J 

/ . \^q ■^q'^p-^p ^pq ^ Sq^p-^p ^pq 

^ ^q '^qSp^pq ^ SqSp^ pq J 
= V r R('=l)r7(^l) + R(al2)r'(a21) _ ( o(a21)^(al2) , t^(c2) n{c2)\\ 

/ J \ pq pg ' pq pq \ pq pq pq pq J J 

p,q 

= trg BC = trg BB""-^ = trg B"^ , (A.3) 

as claimed. 

Appendix B. Details of the Spectral Decomposition 

The matrix K is ordinary Hermitean, although anticommuting variables are present. In 
particular, all inverses of the matrix elements Kn'n exist. The eigenvalues Yn are thus 
uniquely defined. Moreover the diagonalizing matrix = [Vi ■ ■ ■ Vjv] introduced in 
Eq. (|TH|l is ordinary unitary and in SU(A^), and the corresponding eigenvectors Vn are 
orthonormal and have commuting elements only. This might seem to be at odds with 
the form ()25|1 of the eigenvectors Vn. In the second representation of the supervectors 
Wn, the Wnpi are anticommuting and the Wnp2 are commuting, such that all elements 
of the vector Vn are anticommuting, despite the fact that K is an ordinary matrix. To 
clarify this, we use that the Vn form a complete set and expand 

N 

K = E 7nn'K' , (B.l) 
n'=l 

where the coefficients 7„„' are commuting in the first and anticommuting in the 
second representation of the vectors Wn- From the eigenvalue equation we find 

Vi^KVn = YnV^Vn. Inserting the expansion ()B.1|) . we obtain (YJ — Yn)lin = 0, which 
implies that the coefficients satisfy '-fin = Jn^in with new coefficients 7„. Hence we have 

Vn = InVn , (B.2) 

such that the eigenvectors Vn and Vn are proportional to each other. The Vn 
are orthogonal, but they cannot be normalized in the standard way, if the 7„ are 
anticommuting. We emphasize that this causes no problem whatsoever. The coefficients 
can be written as the scalar products 

In = VlVn = V^AL^/^Wn . (B.3) 

If 7„ is anticommuting, it is has to be nilpotent, which means that an integer j exists 
such that 7^ = 0. We notice that 7„ is not nilquadratic, i.e. the number j is here larger 
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than two, because 7„ is a complicated linear combination of nilquadratic anticommuting 
variables. In general, a linear combination of J nilquadratic anticommuting variables is 
nilpotent for every j > J + 1. Moreover, we also deduce from Eq. ()B.2jl 

as the orthogonality relation. 

It is worthwhile to also collect more information about the supermatrix B. 
According to the definition (j^UI) . it is non-Hermitean and satisfies 

= LBL . (B.5) 

One easily sees that 

wl^Lwp = Sp'p (B.6) 

is the corresponding orthonormality relation for the eigenvectors Wp. Being supervectors, 
these eigenvectors can always be properly normalized to unity. The completeness 
relation reads 

2k 

^ WpW^L = hk , (B.7) 
p=i 

where is the 2k x 2k unit matrix. 

We construct a helpful alternative representation of the matrix K. Employing the 
form and the completeness relation ()B.7|1 we work out the expression 

2k 2k 

71 = 1 71=1 

2k 

= ^ AL'/^WnwlLL'/^A^ = ALA^ , (B.8) 

71 = 1 

and by virtue of Eq. (fT^ we arrive at 

2k 

K = Y: V^V^ . (B.9) 

71=1 

This spectral decomposition is somewhat strange, because the eigenvalues do not appear 
explicitly. However, useful results can be deduced from it. In the eigenvalue equation 
KVn = YnVn it givcs together with the orthogonality relation ()B.4|1 

= KlK = 7:77^ • (B.IO) 

Thus, the k eigenvalues 1^+^ = yp2 are products of two nilpotent anticommuting 
variables. 

Furthermore, one readily sees from the decomposition ()B.9|) that all eigenvalues 
Yn, n > 2k which are different from Yp = ypi and l^+fc = yp2 must be zero. As K 
is Hermitean, one can convince oneself in the usual way that the eigenvectors Vn to 
different eigenvalues are orthogonal. Let Vn be an eigenvector to an eigenvalue Yn with 
n > 2k. We immediately conclude from Eq. ()B.9j) that KVn = and hence F„ = 0. 
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Finally, we show in an alternative way that the 2k non-zero eigenvalues y„ of K 
coincide with the eigenvalues of B. We write the eigenvalue equation as Bwn = bnW^ 
and consider the orthogonality relation 

= wliL'/yA^AL'/'wr, = wl{L'/yL-'/^L'/^A^AL'/^Wn 

= wl,LBWn = bnwl,LWn = b„6n'n ■ (B-H) 

For n = n' we conclude 6„ = Yn as claimed, and for n n' we observe that the 
orthogonalities of the eigenvectors Vn and Wn mutually imply each other. 

Appendix C. A Supersymmetric Ingham— Siegel Integral 

For 2k X 2k Hermitean supermatrices a and p, we wish to calculate the integral I{p), 
i.e. the Fourier transform (jK!^ of the superdeterminant. As I{p) is obviously an invariant 
function depending on eigenvalues only, we may replace p with r. Up to a certain 
point, we can apply and slightly extend the methods given in Ref ^H] for the case of 
ordinary matrices. Employing the notation of [Appendix A the matrix a is viewed as 



consisting of the element \ the supervector^i = (o2i , • • • , <^if^\ ■ ■ ■ , '^ki^'') 
with k — 1 commuting and k anticommuting variables, the complex conjugate al and 
the {2k — 1) X [2k — 1) Hermitean supermatrices a containing all other matrix elements. 
Because of 

detga" = detga" [a[f' + al{a-)-^(fi) (C.l) 

the integral over a{i can easily be done with the help of the residue theorem. Some care 
is needed, because the bilinear form al{a^)^^cxi is an undetermined complex number 
due to the presence of the imaginary increments. However, as the variables a are 
only parameters in the ai^^^ integration, we may shift the imaginary increments away, 
assuming that a can be inverted. The unitary supermatrix diagonalizing a can then be 
absorbed into the supervector ai. This makes the bilinear form ai{a)~^ai real, and the 
residue is well determined. The integral over the supervector ai is then simply Gaussian 
and we find 

/(p) ~ e(rii)(zrii)^exp(-£:rii) 

J d[a] exp (itrg ra) detg -^^-^^a' (C.2) 

with f = diag (r2i, . . . , Tfci, 'irfc2, • • • , '^?"fc2)- We perform the calculation up to the 
normalization constant which will be determined later on. It should be noticed that 
the presence of the anticommuting variables leads to some differences as compared to 
the corresponding formula in Ref. ^Hj- We can repeat this step k — 1 further times until 
all variables cr^^^-* and all anticommuting variables c^g^^^ and cr^g^^-* have been integrated 
out. This results in 

k 

Hp) ~ n0(^Pi)(^^i'i)'^exp(-£rpi)J(r2) 
p=i 
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Ar2) 



(C.3) 



The remaining integral J{r2) is over the ordinary k x k Hermitean matrix o"^*^^-*. As 
no anticommuting variables appear in the integrand, the inverse superdeterminant is 
identical to the determinant in the numerator. This determinant does not contain 
singularities anymore and thus we dropped the imaginary increments. Upon introducing 
eigenvalue-angle coordinates for a^'^'^^ and applying the Harish-Chandra-Itzykson-Zuber 
integral (j6T|l for U(A;), we are left with an integral over the eigenvalues s^'^\ p = 1, . . . , /c 
given by 

J(r2) ~ J A,,(s(^2)) exp {itTr2S^'^^) det^'-^s^'^^ . (C.4) 

As the Vandermonde determinant appearing in the remaining integral reads 



{c2)^ 



n 

p<q 



(c2) _ ,{c2) 



det 



q-l 



we can do the eigenvalue integrals and have 



Ar2) 



Afc(r2l 



det 



r. N-k+g-l 



(C.5) 



(C.6) 



p,g=l,...,fc 

The Vandermonde determinant A,t(r2) cancels out. To see this one may use the identity 



(-1) 



N- 



AN-k + q-l)\ ,,,g^-^(5(r,2 



<2-^ 



(C.7) 



{N-iy. 

which is easily derived by multiplying the right hand side with a well-behaved, non 
singular function and integrating by parts. Thus, we have 



det 



g^-^+'?-^(5(rp2) 

a Af-fc+g-l 

orp2 



p,q=l,...,k 



det 



det 



,-,g^-^^(r,2) 



p2 



k—q 
'-p2 



P,q=i,...,k dr^-i 



A.(.2) n '^"'^'^'^ 



p=i 



dr^2-' 



(C.8) 



Collecting everything we arrive at the first of the results ()54p . The normalization 
constant c^k is found by integrating J(p) together with a normalized Gaussian, 



C 



Nk 



{-iy-\N-l)\ 



i2'K 



= /(O) = 2'^('=-i) j d[p\ exp (-trgp2) /(p) 

= j d[a\ exp (-^trg a^^ detg "^(J- = 2-'^('-i) . (C.9) 

When calculating /(O), we used 0(0) = 1/2, as follows from a careful inspection of 
the step leading from Eq. ()C.1|) to Eq. ()C.2|) . The second and the last equality sign in 
Eq. (frT9|l are due to the Efetov-Wegner-Parisi-Sourlas theorem P EH 122 EHl EH , and 
the equality of the two integrals is a direct consequence of the definition 
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Appendix D. Properties of the Fundamental Correlations and the 
Fundamental Kernel 

The fundamental correlations (jU^ are closely related, but not identical to the 
superdeterminant, i.e. to the Fourier backtransform of the function /(p), defined in 
Eq. ((SSI)- This is so, because the Efetov-Wegner-Parisi-Sourlas contributions Pfl EH 
122 123 121] for the function J(p) and for the whole generating function Zk{x + J) are 
different. To gain further insights, we rewrite the fundamental kernel (jTlj) . We introduce 
two 2x2 Hermitean supermatrices a and p with eigenvalues s = diag (spi, isg2) and 
r = diag (ri, zr2), respectively. We then can cast Eq. (jTlj) into the form 

^^'^'^'^Hv, ^^.2) = - , - V . ( / d[p] exp (-^trgpa) I{p) - l) , (D.l) 

7r(Spi - lSq2) \J J 

which is easily verified with the help of formula (|65|). Importantly, we have to subtract 
the Efetov-Wegner-Parisi-Sourlas contribution, i.e. unity in the present case, from the 
integral. The integral over p is now precisely the Fourier backtransform of Eq. (j53p and 
we have 




- isq2 



[D.2) 



which is the result (I73|) . We notice that the fundamental kernel and the fundamental 
correlation function coincide for k = \. For > 1, this is not so and the calculation 
just presented works for the fundamental kernel, but not for the whole fundamental 
correlation function, because we do not know the corresponding Efetov-Wegner-Parisi- 
Sourlas contributions explicitly. 

To calculate the double integral (f7T|) directly, we apply techniques developed in 
Refs. [T^ 123] . We integrate the r2 integral — 1 times by parts and have 



+00 



— 00 

d^~^ exp {-ir2Sq2) 



drn ri — ir2 



(D.3) 



r2=0 



The — 1 fold derivative is computed with Leibnitz' rule according to 



exp {-ir2Sq2) 



' 2 

which yields Eq. ^ 



dro ^ ri - ir2 



N-l 

= (-if(N - 1)! y , (D.4) 

r2=0 n=U \ 
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Appendix E. Rederivation of the GUE Correlation Functions 

Again, we resort to techniques developed in Refs. ^21 12^] • For a Gaussian probability 
density P{H) ~ exp(— trif^), we obviously have 

P^''-^-^\H,,) = ^exp{-Hl) (E.l) 

in Eq. (fTTj) . We plug this into Eq. (fTIj) and use the first of Eqs. (f?^ . The H(^k+q){k+q), Q = 
1, . . . ,k integrals are then recognized as representations of the Hermite polynomials 

on +~ 

H^ix) = — y exp (-e) {z - , (E.2) 



while the Hpp, p = 1, . . . , k integrals yield the generalized Hermite functions 

H^{x) = exp (x^) J exp (-^^ - t2x^) Cd^ . (E.3) 



These functions have already been introduced into the RMT context in Ref. 
They comprise both fundamental solutions of the Hermite differential equation, the 
polynomials Hn{x) and the non-polynomial solutions Hn{x), in the form 

Hn{x) = Hn{x)+tHn{x) , (E.4) 

details can be found in Ref. j2SI- We also use the generalized oscillator wave functions 

1 f 

(pn{x) = , exp Hn{x) = lfn{x) + i^n{x) • (E.5) 



Hk+q){k+q) 



'2"n!v/vr 
Collecting everything, we find 



= exp 
where 



V 1 f 1 

j dHpp—exp(^-H^p^ j dH^k+q){k+q)^exp(^-Hl 

-oo — oo 

C^^'^'^'^\Xp — Hpp, Xq — iH(^k+q){k+q)) 







N-l 

= (pn{Xp)ipn{Xq) (E.7) 

n=0 

is the GUE kernel, including the real part correlations. Hence, we arrive at 



Rk{xi, ...,Xk)= det 



/y>2 ^2 



expi::^^ m^^^\xp,x,) 



■ip,q=l,..;k 



2 ; 

= det [K^'^'^''\xp,Xg) 
= i?f"^)(xi,...,x.) , (E.8) 

which is the correct result. 



p,q=l,...,k 
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